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Abstract. We are interested in the hydrodynamic limit of the Boltzmann 
equation for a gas of fermions in the incompressible Euler regime. We use the 
relative entropy method as improved by Saint-Raymond in the classical case 
|12l . Our result is analogous to the classical case result, but the treatment is 
slightly complicated by the cubic nonlinearity of the collision operator. 

1. Introduction 

The study of quantum gases has been given increasing interest in the hterature 
over the last decade. In particular, quantum kinetic theory is an expanding field 
of research. By quantum gases we mean gases made of quantum particles, that 
are, bosons or fermions. The first ones aim at aggregating together to form the 
so-called " Bose-Einstein condensates" . Conversely, fermions obey Pauli's exclusion 
principle, which prevents any pair of fermions from being in the same state. 

Among the possible models for quantum gases, the Boltzmann-type models first 
proposed by Nordheim in 1928 [TU] then Uehling and Uhlenbeck in 1933 [13] are 
very popular. Although their range of validity is not clear, they seem to capture 
some aspects of the behaviour of bosonic or fermionic particles. Indeed, it has 
been proved, in the bosonic (and space homogeneous) case, that under a threshold 
temperature, a condensate occurs in infinite time |7l [6] . On the contrary, solutions 
of the fermionic Boltzmann equation satisfy a natural L°° bound which reflects 
Pauli's exclusion principle and makes the Cauchy problem easier. There is no need 
for renormalized solutions in this case since one can prove the existence of global 
weak solutions [SJ [5] . 

In the present work, we investigate the incompressible Euler limit of the Boltz- 
mann equation for a gas of fermions (or Boltzmann- Fermi equation). 

1.1. The Boltzmann equation for fermions. It reads in nondimensional form 

(1.1) Maat/ + «.V,/= -^Q(/) 

Kn 

where f{t, x, v) is the density of particles which at time t G M+ are at the point 
X g K'^ of space with velocity w g R'^. The left-hand term expresses the free 
transport of particles in absence of interactions, while the right-hand term takes 
into account the effects of collisions between particles. The Knudsen and Mach 
numbers are positive constants determined by the physical situation of the gas and 
defined by 

mean free path 



Kn 



observation length scale 
bulk velocity 



speed of sound 

The Knudsen number is a measure of the rarefaction of the gas whereas the Mach 
number measures its compressibility. 
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The collision integral Q{f) is given by 
(1.2) 

/ B{v^v,,Lo){f'f:{i-f){i-f,)-fMi-f'){i-f:))dvAuj 

with the usual notations 

= f{t,x,v^), f ^ f{t,x,v'), /*==/(i,x,0 

and where the precollisional velocities (w',w^) are deduced from the postcollisional 
ones by the relations 

v' = V — (v — V^).UJUJ, 
= u* + {v — V*).WCJ, 

Lu being a unit norm vector. These relations express the conservation of momentum 
and kinetic energy at each collision. 

The collision integral (|1.2p differs from the classical one by the presence of cubic 
terms, due to Pauli's exclusion principle, and thus takes into account the quantum 
nature of the gas. Indeed the probability for a particle with velocity v' of taking 
velocity v after a collision is all the more penalized as other particles at the same 
point already have velocity v. This prevents any pair of two fermions from being 
in the same quantum state. As a consequence the solution / satisfy a natural L°° 
bound : if the initial datum /o is chosen such that 

< /o < 1, 

this bound is preserved by the solution / for all times. This property is of great help 
in the study of this equation and implies very different features from the classical 
or the bosonic Boltzmann equation even if they look very similar. 

The function B(z,uj), known as the collision kernel, is measurable, a.e. positive, 
and depends only on \z\ and on the scalar product z ■ uj. It is often assumed to 
satisfy Grad's cutoff assumption: 

(1.3) < B{z,uj) < Cb(1 + \z\f a.e. on x 

(1.4) / B{z,oj)duj > 7^—^ a.e. on 

Js2 Cb 1 + \z\ 

for some constants Cb > and /3 e [0,1]. These assumptions guarantee the 
existence of a solution to equation (|l.ip in the whole domain R'^ or in the torus 
[312]. However, in order to have existence of solutions in more general domains and 
to ensure that they satisfy the local conservation of mass, momentum, and energy, 
as well as the entropy inequality, we will have to make more restrictive assumptions 
[2], as will be explained later. 

The symmetry properties of the collision operator, coming from the fact that the 
transformations (w,w*) i— >■ («',«!) and {v,v') (w*,w^) have unit jacobian, imply 
that, at least formally, 

/ QU)^dv^- j f i? (/'/:(! -/*)-//*(! -/')(!-/:)) 

(1.5) JR3 4 Jr3 Jg2 

and then 

/ {l,v,\v\')Q{f)dv = 0. 

JR3 
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As a consequence, the solution / of equation (jl.ip formally satisfies the conserva- 
tions of mass, momentum and kinetic energy: 

dt J vfdv + ^V,. w ® vfdv^ = 
dt I IH'/d^+^V,.!^! v\v\'fdv^^O. 

Moreover, taking ip ~ log ^ in (|1-5|) leads to the so-called H-theorem, which 
expresses the second principle of thermodynamics: 

dt f s{f){v)dv + ^V,.( f vs{f){v)dv) = -^j^Dif) 

with 



Ma V /b3 / MaKn 



s(/) = /log/ + (l-/)log(l-/) 



and 



Dif) = 7 / Biv~ v^.Lo) (/7:(i - /)(! - /*) - - /')(i - /:)) 
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•^°H//*(i-/')(i-/:)J'^*'"- 

Notice that is non-negative. The minimizers of the entropy, which arc also 

the functions that cancel the collision operator, are given by : 

Proposition 1.1. Assume that g E L^(M^) is such that Q{g) and D{g) are well 
defined and satisfy the bounds < g < 1. Then 

Q{g) = ^ D{g) = ^ g = or g ^ I.^a 

where M is a maxwellian distribution, that is, 

M ~ ae 2b 

with a > 0, b > 0, and for some subset A C M'^. In the first case the distributions are 
called Fermi-Dirac or Planckian distributions, whereas they are called degenerate 
Fermi-Dirac distributions in the second case. 

Proofs and details arc to be found in [5J [7]. The coefhciens of the Planckian 
distribution a,b,u are fully determined by macroscopic parameters of the fluid with 
distribution g, namely the mass pg, bulk velocity Ug and pressure Pg, defined by 



Pg{t,x)= / g{t,x,v)dv, Ug{t,x) = — - — - / vg{t,x,v)dv 

and 

Pg{t,x,v)= I \v ~ Ug{t,x)\'^g{t,x,v)dv. 
This relation will be detailed in the next Section. 
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1.2. The formal hydrodynamic limit. In the fast relaxation regime, that is, 
when the Knudscn number goes to 0, the collision process becomes predominant 
and we expect in view of (jl.ip that a local thermodynamic equilibrium is reached 
almost instantaneously. The distribution / describing the gas is then close to a 
Planckian distribution fully determined by some coefficients a, b, u which in turn 
are caracterized by the hydrodynamic fields pf{t,x), Uf{t,x) and pf{t,x). The 
conservation laws then become up to order 0(Kn) 

MadtPf+\7:x-{pfUf) =0 

MadtipfUf) + ■ (yPfUf ® Uf + \pfl^ = 

Madt{pfu)+pf) + \Jx ■ (^pf\uf\'^Uf + ^PfUf^ = 

which are the compressible Euler equations for perfect gases [I^ stands for the 
three dimensional identity matrix). If now the Mach number also goes to 0, that 
is, in incompressible regime, the first equation becomes V^; • (pfUj) = 0, which is 
nothing but the incompressibility constraint. The other equations of motions are 
obtained by a systematic multiscale expansion, depending on another parameter: 
the Reynolds number, defined by 

Ma 

for perfect gases, which measures the viscosity of the gas. In all the sequel, we are 
interested in the inviscid incompressible regime, so that we consider 

Ma = e, Kn = £«, q> 1 

and investigate the asymptotic e — > 0. In this scaling, Kn = e^^'^ — oo. The 
Boltzmann equation for fermions now writes 

(1.6) edtf + v.\/xf^^Q{f). 

Moreover, the hydrodynamic fields will be assumed to be fiuctuations around a 
global constant equilibrium {po,0,po), so that, denoting the fiuctuations around 
the mass, momentum and pressure by p, u and p, 

pf ^ Pn+ ep, Uf = eu, pf = po + ep- 
We also define the temperature Tf by 

Pf = PfTf 

and its fluctuation T by Tf = Tq + sT where Tq is a constant equilibrium value. 
Plugging these expressions into the hydrodynamic equations (|1.2p we get at leading 
order 

V,u = 0, V,(Top + por) = 

which are known as the incompressibility and the Boussinesq relations; at next 
order it comes 

dtp + u- VxP = 
dtu + u- VxU + Vx{pT) ^ 0. 
The challenge is now to make rigorous this formal limit. 

The rest of the paper is organized as follows: in Section [21 we recall some facts 
about equation (|1.6p . In Section |3] we explain our strategy which is based on the 
modulated entropy and we give our main result (Theorem 13. ip . In Section 0] we 
compute the time derivative of the modulated entropy. In Section [5] we construct 
an approximate solution of (jl.6p the parameters of which satisfy the incompressible 
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Euler equation. We give some useful intermediate results and estimates in Section 
El Finally in Section [7] we end the proof of the main result by controlling some flux 
terms. 



2. Some details on the Boltzmann equation for fermions 

In this Section we explain the relation between the hydrodynamic fields and the 
coefficients a,u,b of the Planckian distribution defined in Proposition p.ip . We 
then state the existence result to be used in the sequel. 

2.1. Relations between the coefHcients and the macroscopic quantities. 

We recall here the definition of the hydrodynamic fields: for a given distribution 
function / the mass pj, the bulk velocity Uf and the pressure pf are defined by 

Pf{t,x)=[ f{t,x,v)dv, uf{t,x)^ — ^ — -/ vf{t,x,v)dv 

Jr3 Pf(t,X) Jg^i 

and 

Pf{t,x,v)= / \v -Uf{t,x)\'^f{t,x,v)dv. 

When the thermodynamic equilibrium is reached, the density function is an equi- 
librium solution 

/= , M = ae ~, 

with a, 6 > 0. It is then clear that u = Uf. Let us define for comodity the functions 

-'R^ 1 + ae-T- 
They are well defined for a e [0, +oo). Then we write 

Pf = 5^/2^0 (a), Pf = 6^/2^2 (a). 
If we define now the internal energy 

we can quote the following result from [3]: 

Proposition 2.1. There exists a positive constant I such that the mapping 

(a, 6) G (0,+oo)2 ^ (P/,e/) G 
where £ = {(p/,e/) s.t. e/ > Ip^j^}, is one-to-one. 

From now on, we will denote by Pq a global equilibrium defined by 

(2.1) Po = T— TT' ^^^0 = e°V"W 

where oq , &o > and oq is such that 

F2(e"«) > 2^^0(6°°)*. 

This ensures that perturbating ao and bo will leave the corresponding perturbed 
mass and internal energy inside £. 
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2.2. Existence theory. Let be a subset of regular enough such that the 
normal is well defined on the boundary. The equation 

(2.2) dtf + v.V,f^Q{f) 
must be supplemented with an initial condition: 

(2.3) /(O, X, v) = fo{x, v) y{x, v)enx W", 

and a boundary condition we choose to be specular reflection for simplicity: 

(2.4) f{t, X, v) f{t, X, R^{v)) V(x, v) ednx such that n{x).v < 

where n(x) is the outer unit vector normal to the boundary and Rxiv) is the specular 
reflection law 

(2.5) Rx{v) =v-2{v.n{x))n{x). 

The choice of the specular reflection as a boundary condition makes all the boundary 
terms vanish in the weak formulation of the equation. It also cancels the Prandtl 
layer along the boundary. 

The following existence result was proved in [51 H] : 

Theorem 2.2. Let fl be either or a regular subset ofMP. Let the collision kernel 
B be such that 

(2.6) < B e L\R^ X S^), 
and let 

(2.7) ,heL°"{nxR^), 0</o<l. 
Then, the problem (j2.2l) - (j2.4[) has a unique solution f satisfying 

f e i°°(R+ X O X R3), < / < 1 a.e. 

Moreover, f is absolutely continuous with respect to t. 

The assumption made on the collision kernel (|2.6p is very strong, much more 
than Grad's cut-off assumption (|1.3p - (|1.4p which would be enough to ensure the 
existence of a solution [21 . However, it is necessary (at the time of the writing) in 
order to prove that the local conservation laws are satisfied by the solution. It is 
the case if we assume that the collision kernel has the symmetry 

(2.8) B{w,Lo)=q{\wl\w.u:\), 

which is physically relevant. The following proposition was proved in [5] : 



Proposition 2.3. Assume that the collision kernel B is as in (|2.6p and 

Assume moreover that the initial datum is as in (j2.7p . and 

(1 + |w|^)/o(a;, v)dxdw < +oo. 

axR3 

Then, the solution f to (l2.2D - (j2.4p satisfies, in the distributional sense, the local 
conservation laws 



dt / /d« + V,. / vfdv] =0 
dt w/du + V^.(/ v®vfAv\=Q 



(2.9) { 

JR3 

dt I Kf /di; + V,. ( / v\v\''fdv] =0 

JR3 \JR3 / 

Finally, under the same assumptions on B, Boltzmann's H-theorem is true: 
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Proposition 2.4. Let the collision kernel B satisfy (|2.6|) and (|2.8|) . Assume that 
the initial datum satisfies (j2.7[) . and 

(|a;p + |wp)/o(x,?;)dxdi; < +00, 

and Zei / be the solution to (|2.2p - (j2.4p . Then, 



s{f){t)dvdx+ / / Z?(/)da;d,s= // ,s(/o). 
axK3 Jo Jn JJnxR^ 

2.3. The linearized collision operator and the fluxes. It is sometimes inter- 
esting to linearize the Fermi-Dirac collision operator Q{f) around a global equilib- 
rium. Let — ^^^^j be a global equilibrium depending on the parameter e, with 

Afg = a^e ^ . Then the linearized eoUision operator writes 



R3 



1 + M,., 1 + A/' 1 + M. 



X (g(l + A/,) + + ^^e.*) - .9'(1 + M^) - g'Al + J) d«*dw. 

The properties of this operator will be detailed and proved in Section [SI Let 
us just mention that if the collision kernel B satisfies Grad's cutoff assumptions 
(|1.3p - (|1.4p then is a Hilbert-Schmidt operator with kernel 

KerLg = span < , , 

\ 1 + 7\4 1 + 7\4 1 + A/, 

This property will be crucial in the proof of our main result. 

Remark 2.5. Let us make here an important remark. The collision kernel B 
cannot satisfy both assumptions (j2.6p (required to prove the existence of solutions) 
and (|1.3p - (|1.4p . We will work with two different collision kernels : B^{z,Ijj) = 
l|z|<^ satisfies (j2.6p and will be used in the definitions of the collision operator 
Q(fe) and the entropy dissipation D(f^). 

However, B = \ will be considered when studying the properties of the linear 
operator. 

In proving the hydrodynamic limit, the collision kernel B^ will naturally appear 
in expressions to control. When the linear operator Lg will be needed, we will write 



Be^ B + [Be - B) 

and use the linear operator properties on the term containing B , while the other 
one will be controlled by hands using the fact that B^ and B are equal on a very 
large ball. 

Let us define the energy and heat fluxes 

= ^ {{v - eu)®"^ -\v- eup/) 

where is the unique number such that ji^^F ^ (Kcr/Zg)-^ and / is the 2x2 
identity matrix. These quantities naturally appear in Section |4l Taking e = 0, we 
can prove that 

e (Ker£o)-^ 

where $ = $0 and 5* = 'I'o- Since £0 is a Hilbert-Sehmidt operator, we can define 
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(2.10) l> = /:o^$ and * = 

$ and Vl/ will appear in the control of the flux terms. 

3. Strategy and main result 

In investigating the ineonipressible Euler limit of the Boltzmann- Fermi equation, 
we use the modulated entropy method as in the classical case [12]. Let us first 
introduce the relative entropy of the solution to (|1.6p with respect to the global 
thermodynamic equilibrium Pq (defined by (j2.ip '): 

+ {l-.fe) log - (1 - A) + 1 - Po) dvdx. 

It is a non-negative Lyapunov functional for the Boltzmann-Fermi equation thanks 
to the H theorem. In the fast relaxation limit, is supposed to be close to a 
thermodynamic equilibrium, we therefore define the modulated entropy by 



Pe 

+ (1 - fe) log 1^ - (1 - /,) + 1 - P,^ dvdx. 

where P^ ~ i'^^m some local equilibrium which approximates the solution /g, 

with = a^e ^ . The relative entropy functional measures this approxima- 
tion since 



H{fe\Pe)> f (^e - dvdx, 



inx 

consequence of the pointwise inequality 

{Vx - \/y)^ < x\og- - X + y \fx,y>0. 

y 

Our strategy consists in studying the time evolution of -ij H{f^\P^). This has the 
good scaling since we want to observe the fluctuations of around an equilibrium. 
To do that we 

• determine the equation satisfied by ■^H{f^\Pg), using the local conservation 
laws satisfied by J^; it contains an acoustic (fast oscillating) term and flux 
terms 

• control the acoustic terms by specifying a^, and 

• control the flux terms thanks to a priori bounds on the linear collision 
operator. 

Our main result is the following: 

Theorem 3.1. Letil be some regular bounded domain o/M"^. Let {fs,in) be a family 
of measurable nonnegative functions over 17 x R"^ satisfying the bounds 

{l + \x\^ + \v\^)f,,,„dvdx<C„ 

la JR3 

< fe,^n < 1, 

and the scaling condition 

(3.1) ^H{f,,„,\Po) < C. 
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Without loss of generality, assume that the fluctuation ge,in defined by /g^, 
Po(l + £ge,tn) converges 

1 + Mo V oo bo 2&0 
We assume moreover that the initial data is well-prepared: 

(3.2) ^H{f,^„^Pe,^n)^0. 

Let be some family of solutions to the scaled Boltzmann- Fermi equation 



(3.3) 



edtfe + v.^,.f, = j^Q{f,) 
fs{0,x,v) = fsAn{x,v) onflxM.^ 

with q > I, and with a maxwellian truncated collision kernel B^(z,lo) — l|j,|<^, 
endowed with the boundary condition 

(3.4) /j(t, a;, v) = f^{t, x, R^{v)) V(x, v) e dVl x such that n{x).v < 

where Rx{v) is the specular reflection law (|2.5p . Assume that we have the following 
control on the tails : 

2 



(3.5) / Po ( „ ) dw < C a.e. 

Jk3 



fe_Pa 

. Po 

and, finally, assume that 

(3.6) $ = and ^ = £q are at most polynomial as \v\ — >■ oo. 

Then the fluctuation (g^) defined by f^ = Po(l + e.9e) converges in i^([0,T]; (1 + 
|fp)Podwdx) weafc 

1 fa u b 2 

where (a, u, 6) is the unique Lipschitz solution to the incompressible Euler equations 
dtu + uS/ xU + V xP = 0, V x-u = 

supplemented with 

9ta + u.Vxfl = 0, 9i6 + S.V:c& = 0, Vx(&oa + to&) = 

on [0,T], and T is the maximal lifespan of the solution. 

The result is not very surprising since it is very similar to the classical case |12j . 
Assumption (j3.6|) was proved in the classical case [8], and it should be possible 
(although very technical) to prove it also in the fermion case. 

Assumption p.2p allows us to consider solutions which are almost instanta- 
neously at thermodynamic equilibrium, thus avoiding the Knudsen layer. To get 
rid of it we need to get a better understanding of the relaxation mecanism. 

As to assumption p.Sp . it is a purely technical control of large velocities which 
allows us to treat flux terms. The exact same assumption is done in the classical 
case and it seems to be a great challenge to avoid it. It has been proved that some 
regular solutions to the classical Boltzmann equation satisfying p.5p exist, but such 
solutions in the fermion case are still missing to the author's knowledge. 

The cubic nonlinearity of the collision operator adds a lot of technical difficulties 
in the proof: a lot more terms, respect to the classical case, appear and must be 
controlled. The L°° bound is of great help to treat these difficulties, but does not 
bring more, so that at the end we do not recover a better result than in the classical 
case. 
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4. The modulated entropy 



The main idea to prove the convergence of the solutions to the scaled Boltzmann- 
Fermi equation toward a solution of the incompressible Euler equation is to study 
the time evolution of the modulated entropy. Let 



1 + A/e 

be a local equilibrium which approximates the solution /g, and take 



(4.1) 



Recall that oq and fog are positive constants which were defined in Subsection 12.11 
ai, u and &i are regular functions of {t,x). 

Proposition 4.1. Any solution to the scaled Boltzmann- Fermi equation p.6p sat- 
isfies the following identity 

D{fe)dxds 



^H{f,\P,){t) = ^ff(^„|Pe,™) ^ 



fe 1, 



J JVl 



Jo 

V — eu \v — EMp 



2K 



.Ae{ai, u, bi)dxdvds 



\ [ ff Afe,*,.V,6id.Tdi;ds f II fei^e ■ '^.u)dxdvds 



2e 
1 

3^ 



-ebi 



ao+EQi- 



bo 



1 + e 



aa+eai — - 



^''1^2 



■dwda; 



ds. 



where the acceleration operator Ae(ai^u,bi) is defined by 
(4.2) 

(dfO + uS/ xO. 
dtu + u.VxU + bo^^^S/^a + ^{t^ - To)\7xb + ^\7x {boa + rob) 
dfb + u.\'xb+ ^Vx-u 

and To is the limit of as e goes to 0. 

Proof. The modulated entropy can be written as 

Hif,\P,) = Hif,\Po)+ I (/,log^+logi±^^d«dx. 

The H-theorem allows to write 



i/(MPe)(<) = i^(/e,m|Pe,,n) 



1 



(4.3) 



e^+' Jo 



fe - log 



D{f,)dxds 



2b, 



dxdv 



ds 



+ f ff log(l + Ms)da:d^; 

Jo LQ' J Ji2xR3 

An integration per part changes the last term: 



log(l + A4)da;du== i 
nxR^ 3 



1 



ds. 



a^e 



■dwdx. 



Computing the time derivative and using the conservation laws (j2.9p leads to: 
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fe 



OxR3 



dt log Oe + W • Va; log Qe 



£ ( 1 1 

+ — (v — eu) ( dtu + u ■ V^w H — Vj;U : (f — eu) H — ^^sVi: loga^ 



\v-eu\^ (dthe , V^6, 1 , , V^6, 



Introducing the rescalcd translated versions of the momentum and energy fluxes 



*e = 1^ [{v - euf - \v - eu\^l) and = {\v 
and replacing and by (|4.ip we obtain: 



dt 



A - log fle 



26, 



dxdz; 



■ {dtai + u ■ V^ai) 



+ —(v ~ eu) ( dfU + u ■ \/xU -i — beV^fli H — TeV^fei 
b^ \ e e 



\v — euP / 2 
: ( at&i + u ■ y^bi + — • u 



1 



dvdx 



2£ JJiix 

Summarizing, the modulated entropy satisfies 



^ Jo J Jn 

-f 



fe { 1, — I , 7^ ) • Ae{ai,u, &i)da::d?;ds 



f2xR3 



b, ' 2be 

fe (be'i/e ■ ^xbi + $e • Vxu) dxdvds 



1 

7o 



dt 



-ebt 



ao+eai- 



OxI 



1 + e 



-dwdx 



where the acceleration operator is defined by (|4.2p . 



□ 



5. Construction of the approximate solutions 

A global equilibrium solution is not expected to be a good approximation of fe 
in the fast relaxation limit since fast oscillations can take place, such as acoustic 
waves. Hence we have to find correctors in order to obtain a refined approximation 
which will lead to the convenient asymptotics. 
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We want to find Ve ~ (af , m^, &f ) such that Ae{a\, u^,bi) in L^, that is, such 
that Ve is an approximate solution of the system 

(5.1) dtV +-WV + B{V,V) =0, 

s 

where 

V = (ai,w,6i), 
WV = (O, (6oai + Toh) , ^\/,,.u 

w.Vxfli 

B(y,V)^\ u.V^u + bohV^ai + TiV^h 
with moreover the constraint 

.-eti i„i2 



and 



Id fP P"'"! aa+eai ^ 

(5-2) Ae = -,:n -TTTZTd-dx ^ 0. 



Remark 5.1. B is a bilinear operator since we can decompose as = tq + eti + 
o(e) and ti zs a linear combination of ai and bi. 

Such a constrution is done by filtering methods, and we refer to |12j for all the 
details. We will just here mention the main points of the proof of 

Theorem 5.2. Let (a^",u™,&™) belong to H''{Q) for some s > |. Then there 
exists some T > and some family {al'^ ,u'^'^ ,bl'^) such that 



sup lim ||(ai'^,M^'^,6^'^)||L^([o,T];H-(n)) < Cj 
(a^'^'™,w^'^'™,6^'^'") ^ (al",u™,&l") in H%n) as e and N ^ +oo, 
A(a^'^, u^-^, bl'^) ^0 m i2([0, T]xn) ase-^O andN ^ +oo, 

and 

, £.JV 

an+£a, — ^ — ^ 

- _| ,2 '° 

1 l^'l — 



1 + e" 



m i:i([0,T]) as £ ^ and N +oo. 
In L'^ equipped with the norm 

liv^ili^ = Ikilli^ + Ml. + ^\%a, + T,Mh, 

the operator W is skew-symmetric. We can thus define the associated semigroup 
W. If we conjugate (|5.ip with yV(|), it comes 

dt{W {^-^ V) + W Q B{V, V) = 0. 

Defining 
it becomes 

(5.3) 9tl/ + W Q j B ( W ( - J j f, W (--^ V \ = 0. 
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We obtain the first order approximation by taking strong limits in the fihered sys- 
tems. However, the error is not expected to converge strongly to due to high 
frequency oscillations. We therefore need to construct a second order approxima- 
tion. We then need a third order approximation in order to have mass conservation 
at sufficient order. 

5.1. Study of W. We expect the solutions of (j5.3p to have a very different be- 
haviour depending on the spectrum of W. Since it is skew-symmetric, its eigenval- 
ues are purely imaginary, and its eigenvectors satisfy 

WVx = iXVx, 

which implies if A 7^ 



a.6a = -^yb 



2to' 

Vj;6a = «^"A 



^A 



Hence, the operator W has the same spectral structure as the laplacian on fl, 
which means that is is diagonalizable on the orthogonal of its kernel. An explicit 
computation shows that the orthogonal projection on ker W is 

|ri|^^ Jiboa + TQb)dx — bga^ 



no(a, u, fe) = a, Pu 

V •^o 

where P is the Leray projection, that is, the orthogonal projection onto the divergence- 
free vectors. In the sequel. Ha is the orthogonal projection onto ker(M^ — XI). 

5.2. Construction. We just sketch here the construction of approximate solutions 
to (j5.3p . and we refer to [H] for all the details. The method is as follows: 

• Decompose the operator W as W — J2xe<y ^^'-^ plug this into (|5.3p . 
It leads to 

• The first order approximation Vo of y is defined as the solution of equation 
(|5.3p when we take into account only the resonant modes, that is, 

dtVo + J2 "a,, i?(HA,^ Vo, Ha,3 Vo) = 

This equation is known to have solutions in L'^^{[0,T^),H''{fl)) provided 
that V™ G H^{^) and s > |. However, Vq is not an approximation of V 
in the sense that AeiVo) converges weakly but not strongly in to 0. We 
therefore have to add correctors. 

• Hence, we construct the second order approximation plugging V = Vq + 
eV\ -I- o(e) into ()5.3p . To avoid the problem of small divisors, we introduce 
the projection onto a finite dimensional subset of C°°(r2) Jat, which is the 
orthogonal projection onto the N first harmonics of W and the N first 
harmonics in ker W . We then denote = JnVq, and we check that 
/ Vodx = J Vo^dx. We define by 

v,^-Jn E ^|%T^^ST^nA.5(HA,yo-HA.,fo^). 

Afci +Afc2 T^A/s 

Then, Vf^ + eV-^ is an approximate solution of (j5.3p strongly in L^. How- 
ever, it does not satisfy / (Vq -I- eV^)dx — > when e — > 0. 
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• We therefore define the third order approximation Vj^ as in |12] by plugging 
Vj^ = Vo^ + eV/^ + e^Vi^ + o{e^) into dO]), and we check that 

- = + eVf + e'^Vi^ is an approximate solution of ([Q]) . 

- jjj V^^dx when e ^ 0. 

The key ingredient for the above computations to work (see [12]) is that, for 
X, fi ^ with A 7^ /i, we have 

UoBinxV,n^v) = 0. 

Indeed, writing Il\V — {a\, u\, b\), we get 

+ {c^ax + c%x)'^xh^ + (c^a^ + c^b^,)Vxbx 
\ ux.Vxbf_,+Uf_,.Vxbx ) 

but = flp = since A, ^ 0, so that 

1 I ° 

noB(nAy,n^V) = -no liA-V:^^^ + u^,^xux + c^bxV^b^, + c%^,Vxbx 

\ ux-S^xbf^ + u^.S/xbx 



P{ux.VxUf^ + u^.VxUx) 
[ux-Vxbf^ + u^.Vxbx)dx) 

Note that / ux-^ xb^dx ~ ^ |^ / bx-b^dx ~ since \^ [i. Moreover, 

{ux-VxU^ + Uf_,.VxUx)i = (Va:(uA.w^)) + u{ {djul^ - diujj) + u'^ {dju\ - diu'^ 

so that 

noS(nAy,n^y) = o. 

In addition, if A 7^ 0, 



^Ada; ^ [ ^x-Uxdx = i-^ [ u.nddx = 0. 



At the end, we get that 



and 



d 

di 



adx ~ 0, 



which is the key to see that ^ j Vj^ dx when e 0. 

The equation for the non-oscillating part can be decoupled from the rest of the 
system, and writes 

9tnof + noS(noi/,noy) = 0, 

which can be rewritten, with the notation HqV = (a, u, b), 

dta + u.\/ xCL = dtb + u.\/xb ^ 0, \7 x{bf)a + tqI)) ^ 0, 

dtu + u.VxU + 'VxP ~ 0, \/x-u = 0. 

This is the incompressible Euler system. Now it remains to show that the approx- 
imate solution we constructed satisfies the constraint (|5.2p . This is the object of 
the following result: 
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Proposition 5.3. With the above contruction oJV^ = (af ,w , 65^), the quantity 

(5.4) = 4^ // ^^\^\'— %^dvdx 

goes to as e ^ for all N > 1. 
Proof. can be rewritten 



wliere 



and 



Using tlie equality 
we compute 

N _ 'J"0 



1 + ae 2 



2e 

We easily prove the two identities 



J (e^^'f i^2(af )a,6f + ei^^f i^o(af )a*af ) dx. 



and 



which imply that 



Using the fact that A{a^ , , b^) in and the conservation of the mass, we 
deduce that 



0. 



6. Useful intermediate results 



□ 



In this section we define and study the linearized collision operator, and then 
give some bounds that will be very useful in the proof of Theorem 13.11 
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6.1. The linearized collision operator. As we look at solutions of (|1.6p which 
are fluctuations around an equilibrium state, that is, 

it makes sense to use the linearized collision operator, which is defined as 

[ [ , , M.^ 1 1 

Leg = / / B[v - v^,uj)- 



S2 JK3 1 + Me,* 1 + 1 + Mi, 

X ((7(1 + Me) + .g*(l + Ale,,) - .g'(l + M'e) - .g*(l + M^^ ,)) dv.duj. 

Let us note 

[ f B{v~V,)MeA^ + Me) , , 

^'■'^ '^^^"^ = i L (l + M.,)(l + M;)(l + Af^.J^^*'^-' 



id 

/e = ttee 21 

£5 with the collision kernel 5 = 1 



and recall that Pe = with Me = a^e ^''^ . Keep in mind that we will use 



Proposition 6.1. Assume that the collision kernel B{z,uj) satisfies Grad's cut- 
off assumptions p.3p - p.4p ; then Ce is a non-negative self-adjoint operator on 
L'^{Medv) with domain 

V{Ce) = {g e L^iMedv) I l^eg G L^iMedv)} = L^il^eMedv). 

and kernel 

K r f 1 ^ 
KerLe ~ span < , , 

^ \l + Me I + Me I + Me 

It can he decomposed as 

Ceg{v) = Ue{v)g{v) + ICeg{v) 
where K-e is a compact integral operator on L^{Medv) and 

< l/_ < Veiv) < + \vf) 

provided that 

(6.2) -^<ae<Co and^<be< Co 

with Co a positive constant independent of e, and v+^v- depend only on Co and 
Cb from dLll) (for B = 1 we have = ^^^^'g', and v+ = 47rCo^^;. 

Remark 6.2. In the framework of Theorem \3.1l assumption (|6.2p is satisfied, since 
we consider for Oe and be small perturbations around the constant values e"" and 
bo. 

Proof Let g, h G L'^{Medv). Then 



f f Me 

J hCegMedv = J hil + Me)CegY:^dv 



I If f , , M. Me, 1 1 

= -// / B{v-v,,u) : '-^ 

477752 ^ 'l + Mel + Me,,l + M'^l + Mi, 

X (g(l + Me) + g,{l + Me,,) - g'{l + M'^ - .g^l + A/;,)) 

X (/i(l + Me) + h,{l + Me.,) - h'{l + Me) - h',{l + M^ J) dwdw.dw 

which easily shows that Ce is a non-negative self-adjoint operator. Moreover, letting 
h = g implies that g is in the nuUspace of Cg if and only if, for almost all (v,v,,ijj) G 
M3 X ]R3 X 5*2, 

5(1 + Me) +g,{l + Ale,,) = .g'(l + A^e) + 9', (1 + K*)- 
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In other words, g(l + Mg) must be a collision invariant, and therefore the nullspace 

of £e is made of all linear combinations of j^f^T' i+m ' i+m (^^^ El)- 
The next step is to split into two parts: 

with defined by (|6.ip and 

M,. 1 1 



^eg{v)= / / B{v - v^.u:) — 



+ M,^, 1 + 1 + Mi, 

X (541 + 14,*) - g'{l + AQ - ffUl + J) dv.du. 

The bounds on ly^ can be proved using Grad's cutoff assumptions (jl.3p - (jl.4|) . They 
are uniform in e thanks to (|6.2p . The operator /C^ can be split into two operators 
in the following way: 

with 



Jr3 Js2 1 + MP + 



£1* 



and 



1 + A/e,* 1 + A/^ 1 + A/; , 
X (^'(l + A//;) + 51(1 + A'/; J) dz;*da;. 

In the classical case, the same decomposition holds, and a clever change of vari- 
ables known as "Carleman's parametrization" allows to show that the following two 
operators are compact on L^(Afedw) (see [TT] for example): 

^ig{v) = I / -B(w — w*,a;)A/e,*g*dw*dcLi 



K.2g{v) = / / B{v - v,,uj)M^^^.{g' + gl)dv^.duj. 
It is easy to see that 

ll^effl|L2(M,du) < ||A.g||L2(A/,dtO 

and 

\\^eg\\L2{M,dv) < \\!^2g\\L2{M,dv) + 11^2 (A/e^) || L2(M^dD) 

from which we deduce that K-^ is a compact operator on L'^{Al^dv). □ 

With these results, we can assert that is coercive, and therefore is a Fredholm 
operator: 

Corollary 6.3. There exists Cc > such that, for each g e T>{C^) n {Ker{Ce))'^ , 

J gCegM.dv > CcWgWh^^^M^dv}- 
If (j6.2l) holds, then Cc depends only on Co and Cb, but not on e. 



Proof. This result is a consequence of Proposition lOl Indeed, the multiplication 
operator g t-^ v^g is self-adjoint on L^(i/£ Afgdu) and has continuous spectrum, 
namely the numerical range of [v^ , -^oo [ where 

V- < ^ hrf Vpiv). 

Then, since JCe is self-adjoint and compact on L^(Afedw), Weyl's theorem ensures 
that the spectrum of Ce is made of \vj , +oo[ and of a sequence of eigenvalues on the 
interval [0, z^~] with ly" as unique possible accumulation point. Consequently, there 
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exists a smallest eigenvalue which is bounded by below by if (|6.2p holds, and 
the following spectral gap inequality is satisfied for each g £ Vi^Ce) n (Ker(£e))^: 

j gC,gM,dv>Cc\\g\\l.i^M^Avy 

Recalling that 

gCegM^dv ~ I v^g^M^dv — / glC^gM^dv 



and using the continuity of /Cg, we get the inequality with the weighted norm as 
stated above. □ 

6.2. Useful bounds. This subsection lists some a priori bounds which are needed 
for the proof of convergence in the hydrodynamic limit. The first estimate comes 
from Young's inequality: 

Proposition 6.4. For z > — 1 we define the function 

h{z) = (1 + z)log(l + z)- z. 

It satisfies 

h{z) > for z > -1 

and 

p\z\<Xh*{p) + jh{z) 
where h* is the Legendre transform of h 

h*{p) = max(pz — h{z)) ~ e^ — p — I. 

2> — 1 

The study of the function h is motivated by the relation 

H{f,\Pe) =11 Peh {^-^^\ d^dx + // (1 - P,)h {^f-jA dvdx, 

Proof. The first property is immediate from the definitions of h. The second one 
comes from Young's inequality 

pz < h*{p) + h{z), 
supplemented with the two inequalities 

h{\z\)<h{z) z>-l 
h*{Xp) < X^h*{p) p>0, A e [0,1]. 

□ 

It will be useful to work in since it is the natural space for the study of C^, 
and then we will use renormalizcd fluctuations instead of the natural one dcflncd 
as fe = Pe{l+ege): 



Proposition 6.5. Let us define the renormalized fluctuations 
and 



e \\l 1- P, 
Then 

(6.3) / / Peg^^dvdx + [[ (1 - PShldvdx < l-H{f,\P,). 



OxR3 
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Moreover, if (j6.2p holds holds, 

(6.4) WsVP'egeh'^idvdx) + \\£k\\L'>-{dvdx) < 4 + y^. 

Proof. Estimate (j6.3p is a direct consequence of the definition of and and the 
following inequality: 

x\og^ - X + y > {y/x ~ y^)'^, Vx,y>0. 

Recalling that < < 1 and that is bounded indipendently of e the second 
point comes by direct inspection. □ 

In the setting, the collision kernel needs also to be renormalized. We repeat 
here the definition of the entropy dissipation: 

Proposition 6.6. Define the renormalized collision kernel by 
with 



(1 + M,)(l + Af,,,)(l + M',)il + M^A' 
Then 

\m\U.,^p^^^^<Dif,) 

where is the collision frequency defined by (|6.ip . 

Note that the scaling condition p.ip and the H theorem imply the following 
bound on the entropy dissipation 



-^11 D{f,)dxds < a 

--^ Jo Jn 



10 Jn 

Proof. This estimate is obtained by the Cauchy-Schwarz inequality: 
with 



We easily see that 

r I 

B^K^dv^duj = 



1 + Me 

Next, using the classical inequality 

{x - y)\og- > AAx ~ ^/yf , Wx,y>0, 

y 

we get the result. □ 

Mixing together the previous estimates, we can prove a relaxation result. Define 
He as the orthogonal projection in I/^(Pjdw) on Ker C^. The following proposition 
will be useful in the control of the flux terms: coupled with the control on the tails 
p.Sp and an interpolation argument, it will allow the control of the third moment 

Of/e. 
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Proposition 6.7. Define and as in Provosition \6.5\ and assume 



Then 



ge- he 9s- he 



1+Me ' 1 + Mr 



Pe 



Proof. Plugging the identities 



into 



gives 



fe = PAl + ege)\ 1 - A = (1 - P.)(l + ehef 



^Jf'J'eA^-fe){l~fe,.) 



wliich leads to the key following identity: 



(6.5) 



+ e^Kl/^(g'eg'e,, - ge9e,* + hek,* - KK,* + 9'ehe - 9eK 
+ g'ehe.* - geh'e,^ + g'e,*he - ge,*K + 9'e,*he,* ~ ge,*h'e,, 
+ e^Kl/^{g'X^,{he + he,.) ~ geCJeAK + K*) 



+ hehe,.{g'e+g'e,:) " h'X,M + 
+ i^*Al/2 {g'X Xk,* - gege,*Kh'e,* 

Raising it to the square and dividing by we get, thanks to 
(6.6) 

1 f 9e- he ge~ K 

< ^/// {sJnfU^ - fe){l - feA - sJfefeA^ - - f^f dvdv.du 

where C depends on Cq. Note that here we used the fact that, for a given integrable 
function u, we have (see [T] for more details) 



u{v')dvdoj < C / u{v)dv 



The next step is to decompose B = Be + {B — Be), so that 



< ^D{fe) + ^ I I l{B~ Be)fefe,*dujdvdv, 
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Recalling that /g can be written = Pi.{l + 2ege + e^gl) and = a^e ^'•^ , 
and thanks to the following inequality 

(6.7) l|„_,^|>^MeM,,, < al/^M^/^Mli''e'^^ , 

we get 

C 



£2 



\v—V:^ I > -\ 

/e/e,*da;dudw* 



+ '^^'^ jjj hv-v,\>^PePe.*9%,*dvdv,+Ce'\Ml2^p^^,y 

Using Proposition 16.41 we obtain for 5 > 

<S I I P,.M I ^^-A] dvdv, + ff ,^ 1 P,P,.,/i* ( -X 1 dvdv,. 



We choose 6 = 20Coe^ : 

.PePe,*gldvdvt 

1/2 

jj p,..Peh(k-Jiyi^Av^^ 

We now apply the following variant of Young's inequality 



0(1). 



< I ^ ) +/3r, y > 0,/3 > 0,a e (0,1) 



with /3 = 1, a = i and y = (j^ JJ P,,,P,h [^^) dvdv^Y : 
l|„_„,|>4.P£Pe,*.gedwdt;: 



1 /"/"^ r.uffe-Pe 



2 



< ce'^ (^^ yy (^'^V^ j "^'"^"v + 

Since Sy/ P^^^ge,* is bounded, the same computation works for the second term 
in the right hand side of (|6.8p . so that we finally get 

^JJJ B [^rJU^ - .fe){l - fe,.) - ^.fefeAl - - fi^.)f du^dvdv.. 

since -pD{f^) = o(l). Pluging this inequality into (j6.6p and using the coercivity of 
the linearized operator fCorollarv l6.3p gives the result. □ 

7. End of the proof 



We now have all the ingredients to prove Theorem 13. II We begin by controlling 
the flux terms, and then end the proof of the theorem. 
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7.1. Control of the flux terms. Our goal is to estimate the terms 
1 /■* 



JJn 



and 



4 / // fei'S>e-^xU)dxdvds 

£ JO JJnxs.^ 



with respect to the modulated entropy and the entropy dissipation. The main 
difficulty is that Vl/e = 0(|wp) whereas the modulated entropy allows us to control 
the powers of v only up to 2, via Young's inequality. We therefore try to gain as 
much integr ability as possible. Here the relaxation estimate (|6.3p plays an central 

role. Since it gives a control on fq^^g^, we use a new decomposition of 

(7.1) ^ = P^+^'Pe(l^](9s + h ■ 2 



1 + M, V e 



and we obtain the following estimate: 



Proposition 7.1. Assume that (f>e and ^'^ have at most polynomial growth and 
that (|6.2[) holds. Then the flux terms are bounded by the following quantities: 

4/ // febe'^e-'^xbidxdvds + ^ f ff • V^w)da;di;ds 

^ Jo jJnxs.^ ^ JO jJnxs.^ 

H{MP,){s) + j^^ D{f,){s)dx^ ds + o(l). 



c /■* 

-~/ \\Dx{u,bi){s)\\L2^L^f^^-i) 
^ Jo 

Proof. The proof for the momentum flux (involving $j) is identical to the one with 
the energy flux (involving ^e), hence we will focus on the first one. Writing as 
in (j7.ip . we use a first decomposition of the flux terms: 

(7.2) 1 / <^,f,dv =- I <!>,Pel^^dv + I <^,P,^^f^{g, + h,)dv 

e J I + IVL J 1 + ' 

since /j^g P^'^^dv — 0. To treat the first term in the right hand side, we use the 
identity (|6.5p . wich implies 

(7.3) 

1 ^ ,..5^^^^^^ / M^^^i^^^e^^^ 



^2 (y^/:/:,,(l~/e)(l-/e,*)- V^/e/e,*(l-/n(l-/L)) dz;d«,dL. 
+ g'ehe,* - 9eK,* " 9e,*K +5e,*^e,* - ffcJl^^*) dwdw^dw 

/Jl B^,A,[g'J'^^,ik + K.) - ge9eAK + K*) 

+ Khe.*{g'e+ g'e,*) " ^i^^k* (Se + ffe,* )) dudw* dcj 

B^^Keig'^g'^^Shehe,* - .gege,*/!^/!^^ Jdwdu^do;. 



The first term can be estimated by the entropy dissipation as showed by the next 
lemma, while the other terms are controlled straightforwardly by the L^{Mi,dv) 
norm of g^ , and hence by the modulated entropy. 
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Lemma 7.2. Assume that has at most a polynomial growth and that (|6.2p holds. 
Then the first term in the decomposition ()7.3[) is estimated by: 

<C ^ ^ D{f,)dx + e^^^"^ 
Proof. Let 

We decompose 



(7.4) 



y S,l>,Ay2G,dwdi;*dw + J {B - B,)^,Ay^G,dvdv^duj. 



The first term is dealt witli tlianks to tiie Caucliy-Scliwarz inequality and Proposi- 
tion dlH): 

J B,^,Ay^G,dvdiuduj <c(yj B{^,fA,dvAvdu}j ' D{J,fl^. 

By Corollarv l6.3[ the coercivity inequality implies, since $e e (Ker£^)-L, that 

j B{^,YAAvAvduj<C J ^^^^{^,)Mv, 

with C > independent of e. 

The second term in the right hand side of (|7.4p . containing the high velocities, 
is handled thanks to the inequality 

using that ^^(l + Me) is at most polynomial: 

JiB- + AQAl^^Gdvdv^duj 



< Ce 32S7pr^ 

□ 



The next terms in ()7.3p are dealt with using Cauchy-Schwarz inequality and the 
a priori bound (|6.4p . and therefore are easily bounded by C'(||5e|||2[p^(j„ + ||/ie|||2)- 
Coming back to (|7.2p . the second term on the right hand side is splitted as follows, 
in order to use the relaxation estimate: 



l + Me J \1 + Me 1 + Me 



(rnfe) 



The last term is easily handled: 

9e - K 
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while we need some control on the high velocities for the first one. From the 
hypothesis p.Sp . we deduce that 

P, (i^) ''dv<Cpa.e. 

J 

for all p < 1 and uniformly in e (which depends on p), since the moments of 
differ from that of Pq only by terms of order e. This estimate leads to some control 
on the rescaled fluctuations: 



e\ge\ + e\he\ = 0{l)L^^Liv(p^Av) 



so that 
(7.5) 



l + Me 



9e 



1 + M, 



We easily find the following interpolation inequality for all functions (j) £ L'^ O L^^ 



which implies, together with (|7.5p and Proposition 16.71 



96- he 5£ - he 



l + Me 



1 + Me 



T2 I 



fe-Pe 



Pe 



dv 



1/2 



= Oy\\ge\\L^P,dv) + Wheh^ 

Choosing p > 5 implies jf+T > 2, and Holder's inequality gives 

il0203-Pedw < ||(?il||L2(P.)||'/'2||^^^^ ^||(/'3|!_|£_(p^) 



0(1). 



such that taking ipi = ije + he, 02 — i'^m' 
final estimate: 



n 



and 03 = $e leads to the 



/ 



1 + Me 



l+Me 



{ge + he)dv = 0[\\ge\\i2(p^^,) + \\he\\ 



Peh 



dv ] +0(1) 

□ 



7.2. Proof of convergence. In view of the preceding results, we are now able 
to prove Theorem 13.11 Proposition 14.11 together with Proposition 17.11 imply the 
following: 



1 



^H{fe\Pe){t) 



1 



1 



£94 



Jn 



D{fe)dxds 



< ^H{,feMPe,.n)-- 



1 

3^ 
C 



/e l, 



J JUx 



V — eu \v — eu\' 



OxR3 



bo 



ao+eai- 



1 + e""^"^"^ 255^ 



2be 



dwdx 



.Ae{ai,u, bi)dxdvds 



ds 



^ \\D,iu,bi){s)U2nL^^n) (^H{fe\Pe)is) + J^D{fe)is)dx^ ds + o{l). 
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We plug the functions (a^'^, w^'^, 6^'^) constructed in Theoreni l5.2l in the approx- 
imate solution P^. It follows by Gronwall's lemma that 

li/(MPf )(t) < li/(^„|n^„)exp ||i?.(u^^^,&I'^)(.s)||L.nL=^(0)) ds 



exp 

1 

e 



fe I 



e N I ^ N\2 



-(- 



2b, 



-dvdx + o(l) 



1 + e °^ 1 ^ 

The last term go to as e — > and iV thanks to Theorem 15. 2[ using the fact 
that 

/ Aeia{^,u'-^^,b{''')dx = 0. 
Jn 

The fist term can be decomposed as 



+ 



P 

If r, M (a'j",ti'",&5") 
Ue,in - -r(a5",«»",b5") j iOg 



+ (-F'(a5",«»",b5") - fe,m) log ^ pN ' ^ 



dwdx 



which tends to as e — > and oo since ■^H{fe.in\P{a\" ,ui",b{")) ~^ by 

hypothesis and (a^'^'™, w'^'^-™, fe^'^'™) ^ (ai",u",61") by Theorem[ 
This leads to 

^H{fe\P,^) ^ in L^^([0; T)) as e -> and ^ cxd, 



which is a stronger result respect to the one stated above. 
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